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' We propose a method to evaluate parameters in the Hamiltonian of the Ising chain under 

f"*^ ' site-dependent transverse fields, with a proviso that we can control and measure one of 

the edge spins only. We evaluate the eigenvalues of the Hamiltonian and the time-evoultion 
operator exactly for a 3-spin chain, from which we obtain the expectation values of o~ x of the 
first spin. The parameters are found from the peak positions of the Fourier transform of the 
expectation value. There are four assumptions in our method, which are mild enough to be 
satisfied in many physical systems. 
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1. Introduction 

Quantum information processing has been studied for a long time since Feynman's pi- 
oneering work. 1-3 Many researchers are working toward physical realizations of a practical 



(N 

o 



quantum computer. We have to overcome several obstacles, however, to physically realize a 
working quantum information processor. One of the main issues is to synthesize a high-quality 
qubit with a long coherence time. If we prepare a set of high-quality qubits, we have to identify 
the interaction strengths among qubits with a high degree of accuracy. Without knowing the 



inter-qubit interactions, it is impossible to control the qubit state at our disposal. In some 
cases, the majority of the qubits must be isolated from the environment to suppress deco- 
herence, in which case we have an access to a single qubit only. Then we need to find the 
parameters in the total system, such as the coupling constants, through the single qubit. 

There are existing works on evaluation of spin-spin interactions with such a restricted 
access. 4-10 Burgarth, Maruyama and Nori, 7 ~ 9 studied the Hamiltonian determination in the 
Heisenberg and the XXZ models with N spins. They took advantage of the fact that the excited 
states in these systems are classified according to the eigenvalues of a conserved quantity and 
each excited state belongs to an iV-dimensional subspace of the total Hilbert space. As a 
result, it is sufficient to work with this A^-dimensional subspace in evaluating the parameters 
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in the Hamiltonian. 

The purpose of this paper is to extend their analysis to more challenging cases in which 
the whole 2 N basis vectors are required to construct an eigenvector. Our particular example 
studied here is the Ising chain in site-dependent transverse fields. We restrict ourselves mostly 
within a three-spin chain, the shortest nontrivial chain, to make our analysis concrete. It also 
turns out that the result is extremely lengthy for a longer spin chain so that it is impossible 
to write down the results in compact forms. 

We introduce our model Hamiltonian and quantum dynamics in the next section. Section 
3 is devoted to Hamiltonian determination. Discussion and conclusion are given in Section 4. 
Details of the calculation are summarized in Appendix. 

2. Hamiltonian and Time-Evolution 

We investigate, in this paper, how to determine the spin-spin interaction strengths and 
the external magnetic field strengths in a 3-spin Ising chain with site-dependent transverse 
fields. The Hamiltonian is given by 

2 3 

n(hi, Ji) = ^ Ji*fa! +1 - , (i) 
i=i i=i 

where of and a\ denote the x-component and the z-component of the Pauli spin matrices at 
the i-th site, respectively. It is assumed, without loss of generality, that the magnetic field is 
applied along the x-axes of spins 2 and 3 by employing local gauge transformations. This is 
due to the fact that the dynamics of of is independent of the direction of the magnetic fields 
h,2 and hs in the xy-plane at the sites 2 and 3, respectively. This can be explicitly seen by the 
identity 

V^afV = of, (2) 

where V = e _ *( a2Cr 2+ a 3°-f) - g a transformation which rotates the spins 2 and 3 in the xy-plane 
by angles ct2 and as, respectively. We take hi and /13 to be positive by making use of this 
freedom from now on. 

Our purpose is to evaluate the parameters in the Hamiltonian by controlling and measuring 
the first spin only. We assume the following four mild conditions to this end: 

• It is possible to prepare the fully polarized state | ttt) as an initial state. 

• We can control the first spin only. We may prepare the states such as (a| t) + P\ TT) 
by making use of this assumption. 

• Once the initial state is set up, we can control the magnetic field h\ at the first site only. 
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Other fields fi2 and /13 are fixed but unknown to us. 
• We can measure the spin dynamics of af only. 
Our proposed scheme is summarized as follows: 
Step 1 We prepare an initial state. 

Step 2 The dynamics of the x-component of the first spin (of(t)) is measured. 
Step 3 The Fourier transform of (cj) of (of(t)) is evaluated. 

Step 4 We determine \Ji\ and \Ji\ from the peak positions of of(u>) as functions of 

hi. 

Step 5 There are four combinations of the signs of J\ and J2. We evaluate the short time 
behavior of of(t) numerically for each of the four cases and compare the results with 
that obtained experimentally to determine the signs. (The dynamics is independent of 
the signs of hi and hs as remarked previously.) 

Let {ej}i<j<8 be the eigenenergies of the Hamiltonian, whose explicit forms in terms of the 
parameters {hi, J{\ are given by Eq. (A-l) in Appendix. The eigenenergies of this Hamiltonian 
satisfy the following symmetry relations: 



(3) 



£5 = 


-C4 


= £4 


£6 = 




= £3 


£7 = 




= e 2 


. es = 


-ei 


= e l 



ei + e 4 = e 2 + e 3 = -(ej + e 3 ) = -(e 5 + e 3 ) 

ei - e 3 = e 2 - e 4 = -{e\ - eg) = -(e 5 - e a ) (4) 
ei - e2 = €3 — €4 = — (ej — e§) = —(eg — ej) 
Let us choose the fully polarized state | ttt) as an initial state. The state of the chain at 
time t > is \ij){t)) = U(t)\ ttt)) where U{t) is the time evolution operator 



U(t) = exp(-mt) = ^e-^Vj. 



Here 



v* = n 



k=i 



n-e k i 



(5) 



(6) 
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is the projection operator to the eigenspace with the eigenenergy €j. The real-time dynamics 
{af (t)) of the first spin is calculated with respect to \ip(t)} as 

4 

(af(t)} = (il>(t)\<jf\i/>(t)) = C+Y^ [ A mn cos (e m + e„) t + B mn cos (e m - e„) t] , (7) 

m=l 

where the coefficients C, A mn and B mn are functions of {Jj} and {/ij}. For example, C is 
explicitly written as 

8 

c = (m lE^^tm, (8) 

i=i 

while ^21 and .B41 are 

a 2 i = (ttt I (Piof r 7 + r 2 (r x 1 rs + Vi^Ti + r 8 afr 2 ) | ttt) , (9) 

and 

B41 = (ttt I (^V?^ + Pi + P 8 0?P 5 + W^s) I ttt) , (10) 

respectively. 

It is found from Eq. (7) that the peak positions of the Fourier transform af (oj) are poten- 
tially at e m ± e n for all combinations of m and n (1 < m, n < 4, m / n). It should be noted, 
however, that some combinations of m and n have vanishing amplitudes A^n, B^ n , for which 
case the corresponding peaks do not exist. In the present case, it can be shown exactly that 

A 41 = A 32 = B 21 = B 31 = B 42 = B 43 = °- U 1 ) 

These amplitudes satisfy the sum rule 

C + Y, ( A ™ + B mn ) = (ttt I ( E V 3 J a l ( E V * ) I ^t) 

m=l \i=l / \k=l / 

m>n 

= (ttt I ttt) = 0, (12) 
where use has been made of the completeness relation 

8 

i=i 

3. Hamiltonian Determination 

Suppose we conducted an experiment and measured (af(t)), or equivalently af (cu). Now 
our task is to identify which peak corresponds to which e m ± e n in order to evaluate the 
parameters {</«} and {hi}. Figure 1 shows e m ± e n with non- vanishing A mn and B mn as 
functions of the transverse field hi, which is assumed to be controllable. Other parameters 
are arbitrarily chosen as h 2 = 6, I13 = 7, J\ = 4 and J 2 = 5. The red solid curves and blue 
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hi 

Fig. 1. (Color online) Red solid and blue broken curves denote e m + e n and e m — e n , respectively, as 
functions of the transverse field hi with h,2 = 6, = 7, J\ = 4 and J2 = 5. We have plotted the 
first quadrant part of the graph for simplicity The peaks in the other quadrants are obtained by 
mirror reflections of this diagram with respect to the horizontal and/or the vertical axes. The red 
(blue) dot at hi — shows the nonvainshing amplitude l (B^ j). 

broken curves in Fig. 1 depict e m + e n and e m — e n , respectively. It should be clear from this 
diagram that the curves cross intricately for intermediate hi, while the behaviors of the curves 
for small h\ and large h\ are regular. Then identification of the curves and e m ± e n in these 
regions is easy. It should be noted here that the case h\ = is singular: af (t) remains zero 
for any t > for hi = and IV'(O)) = | ttt)- This case will be treated separately below. 

In principle, we may also use the heights (amplitudes) A mn and B mn of theses peaks to 
determine the parameters. In practice, however, it is rather difficult to use the data for this 
purpose due to the following reason. First, the analytical expressions for A mn and B mn are 
extremely lengthy and we need to employ extensive numerical optimization to determine the 
parameters which fit with the data. Second, the peak has finite width in actual experiment 
due to various reasons, such as relaxation. We must integrate the peak profile with respect to 
u to find the amplitude, which always suffers from error. In addition, reading off the frequency 
can be done with much better precision than measuring the amplitude. Figure 2 shows non- 
vanishing amplitudes as functions of h\ for /12 = 6, h% = 7, J\ = 4 and J2 = 5. Peak amplitudes 
C, A21 and A43 are odd functions of a; while others are apparently not. Note, however, that the 
positive-/ii (negative-/ii) part of A31 and the negative-/ii (positive-/ii) part of A42 combined 
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Fig. 2. (Color online) Peak amplitudes (a) C, (b) A 21 , (c) A 31 , (d) A 42 , (e) A 43 , (f) i? 4 i and (g) B 32 
as functions of /ii. Other parameters are h 2 — 6, h 3 = 7, Ji = 4 and J2 = 5. All other amplitudes 
vanish identically. 



together define an odd function. The loss of antisymmetry in these amplitudes is due to our 
choice of the energy eigenvalues (A-l). If we demand that the eigenvalues be defined uniquely 
by Eq. (A-l) for all parameters, we must give up the antisymmetry of the amplitudes A mn 
and B mn . This also happens to the pair S41 and -B32. Figure 2 demonstrates that the sum 
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rule (12) holds for any h\. 

Let us first consider the limit of a strong magnetic field hi, that is, the parameters satisfy 
\hi\ 3> h2, hs, |Ji|,|J2|. The non-vanishing Fourier peak positions e m + e n (1 < m, n < 4) 
diverge linearly as hi — > oo. In contrast, the peak positions e m — e n converge to constant 
values. We obtain the asymptotic behavior of e m ± e„ for \h\\ » h2, h 3 , |Ji|, from the 
explicit forms of e m given in Eq. (A-l) as 



ei + e 2 - 


■> 2/ii + 


\l (/i 2 + /i 3 ) 2 + A + A 


J(h2 


- h) 2 + Jf 


ei + e 3 - 


2/ii + 


^2 + /i3) 2 + J|-) 


J(h2 


- /is) 2 + A 


£2 + £4 - 


2/ii - 


^2 + / i 3) 2 + «/ 2 2 -) 


J(h 2 


- h 3 ) 2 + Jf 



£ 3 + e 4 ->■ 2/ii - 



+ 
+ 
+ 
+ 



and 



ei — €4 — > 



£2 



2^2 + / i3 ) 2 + J| + °(^)' 
- e 3 -> 2^ (fc - /i 3 ) 2 + Jf + O (J^ . 



From these results, we find the ordering of the peak positions for a large hi as 

£i + £2 > £i + £ 3 > £2 + £4 > £ 3 + £4 

and 

£1 — £4 > £2 — £ 3 - 



(13) 
(14) 
(15) 
(16) 

(17) 
(18) 

(19) 

(20) 



Once the peak positions are identified with particular e m ± e n , we may easily find 
\J (/12 + hs) 2 + Jf and \J (J12 — h 3 ) 2 + J 2 numerically. Clearly these two numbers are not suf- 
ficient to determine the four parameters in the Hamiltonian and we need to examine the peaks 
at other values of hi, such as hi = 0, to acquire sufficient data to completely determine the 
parameters. 

We note en passant that the asymptotic behavior of the peaks does depend on J2 but not 
on Ji. This counter-intuitive character is understood if we study the perturbation theory at 
hi — > 00. Let us separate the Hamiltonian as 



T~L — Ho + Hi, 
Ho = —hio\, 



(21) 
(22) 
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2 3 
i=l 1=2 

The eigenvectors of the unperturbed Hamiltonian %q are ^(|t) ± 1 4)) I s 2 S3}, where S2 ; S3 G 
{t,4}- The first order perturbation of the eigenvalues is independent of J\, which is why 
the asymptotic eigenenergies are independent of J\. This behavior is also understood from 
the spin dynamics. When the magnetic field h\ is much larger than the other parameters, 
the characteristic time r of the motion of the first spin is very short compared to other 
characteristic times. Then the effect of the spin-spin interaction J\ is averaged to vanish over 
the time interval t>T and the J\ interaction is ignored in this limit. Clearly this remains 
true for any number of spins in the chain. 

Next we consider the case when h\ = to acquire sufficient data to fix all the parameters. 
The eigenenergies in this case are ±ei, ±en, each of which is doubly degenerate. Here we take 
ei > en > 0. The explicit forms of eiji are given by Eq. (A- 10) in Appendix. It is easy to show 
that 

{m\u{tf<7iu(t)\m) = Q 

for hi = 0. This is an artifact of the initial state chosen above. Note that af does not commute 
with the Hamiltonian and a proper initial condition introduces a nontrivial expectation value 
of the observable. Let us take an initial state |+ ft); where 

l+) = ^(lt) + U». 

Then the real-time dynamics (af (t)} is calculated as 

(<7f(t)) = C + A° U:l cos (ei + e„) t + cos (ei - e„) t, (24) 

where the coefficients C, Af ll and B® 11 are functions of {Jj} and {hi}, whose explicit forms 
are given in Eq. (A- 15). 

The peak positions with nonvanishing amplitudes at h± = are 

w = ei ± e„ = ^So + 2t/Q ) ±^S - 2^C~ Q , (25) 

and u = 0, where 

S = hi + hi + J\ + J|, C = h\h\ + h\ J\ + jfJl (26) 

Since peak positions ei ± en produce only two numbers y So ± which are essentially 

the same as and Co, we should combine the data acquired for the case of a large hi to 
determine all the parameters /12, /13, \Ji\ and | J2I- Once all the non- vanishing peak positions 
are identified, the measurement of the peaks of af (u>) for additional finite h\ helps us determine 
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the parameters with higher precision. Figure 2 shows that all the peak amplitudes decay as 
hi — > oo. It may happen that the peaks are not observable at large enough h\ since their 
amplitudes vanish in this region. If this is the case, we need to identify peaks at h\ = first 
and then sweep h\ to acquire sufficient data to determine the parameters. 

For the parameter choice, h\ = 0, h 2 = 6, h% = 7, J\ =4 and J 2 = 5, we obtain the 
amplitudes numerically as 

C = 0.5984, Ajj j = 0.05525, = 0.3464. (27) 

They do not agree with the h\ — > limit in Fig. 2 since their initial states are different. These 
amplitudes satisfy a sum rule different from Eq. (12), 

C + A° U j + .Bjj j = (+ tt \of |+ tt) = 1- (28) 

Now that h 2 , /13, I J\\ and | J2I are determined, our next task is to find the signs of J\ and 
J 2 . Figure 3 shows the dynamics of af for | Ji| = 4, | J2 1 = 5, h\ = 1, /12 = 6, and = 7 for the 
four combinations of the signs of J\ and J2. The signs are (sign( Ji), sign( J2)) = (+, +), (— , +), 
(+, — ), and (— , — ) from top to bottom. We immediately note that af(t) is positive (negative) 
for small time limit if sign(Ji) is negative (positive). This is because the molecular field at 
the first site is approximately (— h±,0, J±) for small t. Thus sign(Ji) is fixed by inspecting the 
short time dynamics of (af(t)). Unfortunately, sign(J2) cannot be determined in this way. 
We also note that the dynamics (af (t)) changes the sign under the transformation (Ji, J2) — > 
(— J\, —J2) as seen from Fig. 3. This is due to the following identities, 

3 2 3 2 

X 2 A~ E hot + E J i° Z i°t + i)X{. = - E - E J ^l + i (29) 



i=l i=l i=l i=l 



and 



= (ttt \xl w e in(h >> Ji)t X 2!W afxl 7T e- in{h *' J > )t X 2 , 7r \ ttt) 

where X 2 ^ = exp(— iaf tt/2) is a unitary matrix corresponding to a rotation of the spin 2 
around the x-axis by an angle tt. 

4. Summary and Discussion 

In this paper, we propose a method to evaluate the spin-spin interaction strengths and 
the site-dependent magnetic field strengths only through one of the edge spins. We consider 
a three-spin Ising model with site-dependent transverse fields hi. Since the Ising model with 
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Fig. 3. Dynamics of erf (t) for (sign(Ji), sign( J2)) = (+,+), (—,+), (+,—), and (—,—), where the 
Hamiltonian parameters are hi = 1, /12 = 6, /13 = 7, |Ji| = 4 and \J%\ = 5. 

transverse fields has no good quantum number such as the total S z , which the Heisenberg 
and the XXZ models have, it is considerably more difficult to evaluate the parameters in the 
Hamiltonian in this case. 

First a fully polarized state | ttt) (when hi 7^ 0) or ^=(| t) + | tt) (when hi = 0) is 
prepared and the time-dependence of the x-component of the first spin is measured. The peak 
positions of the Fourier transform af (oj) are at {e m ±e n }, where e n is the first four eigenvalues 
of the Hamiltonian. The identification of each peak position with a particular e m ± e n is easy 
for hi = and hi — > 00, while it is difficult for a finite hi in general. Of course, it is possible 
to identify the peaks for a finite hi if hi is swept from to a large value continuously. 

We determine {hi} and {| J«|} completely from the peak positions of af(u). It was shown 
that some possible peaks are not observable since they have vanishing amplitudes. Some 
combinations of the parameters are read off from the peak positions at large hi and hi = 0. The 
combined data is sufficient to fix all the parameters completely. Additional data obtained for 
finite hi improves accuracy of the estimated parameters. The signs of the spin-spin interactions 
are found by comparing the experimentally obtained spin dynamics with numerical simulation. 

In our proposed method, there are only four mild conditions which are satisfied in many 
experimental systems such as a liquid state NMR quantum computer, which is abbreviated 
as "an NMR quantum computer" hereafter. In fact, the Ising model with site-dependent 



10/16 



J. Phys. Soc. Jpn. Full Paper 

transverse field is a typical system in an NMR quantum computer. It is possible to control the 
magnetic field strength hi of each spin independently by making coordinate transformation 
to a rotating frame of each spin rotating with the respective Lamour frequency, It should be 
noted, however, that there is a subtle problem in an NMR quantum computer associated with 
the initial state. The initial pure state is prepared as a pseudo-pure state, which is produced 
by the time-average or the space-average method, for example. 12-16 To produce a pseudo- 
pure state, however, we need to know the Hamiltonian in advance, which is "the chicken or 
the egg" causality dilemma. The best we can do is to prepare a pseudo-pure state with a 
molecule with a known Hamiltonian and then pretend as if we do not know the Hamiltonian 
and demonstrate our method. We believe our scheme can be experimentally demonstrated 
also with other physical systems with a pure state. 

The authors are grateful to Daniel Burgarth and Koji Maruyama for valuable comments 
and discussions. This work is partially supported by 'Open Research Center' Project for Pri- 
vate Universities: matching fund subsidy from MEXT, Japan. The authors thank the Yukawa 
Institute for Theoretical Physics at Kyoto University. Discussions during the YITP work- 
shop YITP-W-10-14 on "Duality and Scale in Quantum-Theoretical Sciences" were useful to 
complete this work. S.T. is partly supported by Grant-in-Aid for Young Scientists Start-up 
(21840021) from the JSPS and by MEXT Grant-in-Aid for Scientific Research (B) (22340111). 
The computation in the present work was performed on computers at the Supercomputer Cen- 
ter, Institute for Solid State Physics, University of Tokyo. 

Appendix: Derivation of Physical Quantities 

A.l Eigenenergies 



After straightforward but tedious calculation, eigenenergies of the Hamiltonian (1) are 
obtained explicitly as: 

^oA 1 / 2 
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Vdet^ = /if + hi + - J 2 + J. 



2\2 



-2/i 2 (/i 2 , + h\- Jf + J|) + 2hj (-hi + J 2 + Jf ) . 

Figure A-l shows the eigenenergies as functions of the transverse field h\ with J\ = 4, 
J2 = 5, /12 = 6, and /13 = 7. All eigeneneries are doubly degenerate when h\ = 0. 

The eigenenergies €4 and 65 cross when £4 = £5 = 0, which seems to be in contradiction 
with the no-level crossing theorem. 1 Nonetheless this is allowed since the eigenvectors cor- 
responding to these degenerate eigenenergies are symmetric and antisymmetric states with 
respect to the spin indices. The symmetric state and antisymmetric states are characterized 
by introducing the parity operator: 18 



X = af <g> of <g> erf. 
The Hamiltonian and the parity operator commute: 

[H,X] =0. 

Symmetric and antisymmetric functions are defined as follows: 

\^ s ) = Y^a s (\s) + X\s}), 

W 



(A-3) 
(A-4) 

(A-5) 
(A-6) 



|*as} = 5]'a a (|s> - > 
W 

where \s) = \s1S2S3} denotes a spin configuration, and Yl'{ s } denotes summation over all 
spin configurations with si = +1. Observe that Af|^ s ) = \ty s ) while Af|^ as ) = — |^ as ). The 
eigenstates with the eigenenergies £4 and £5 belong to different symmetry classes and their 
eigenenergies may cross as h\ is changed. 



A. 2 hi = Case 

It was shown in the main text that the measurement of the peak positions of af (uj) at 
a large hi was not sufficient to determine the parameters in the Hamiltonian. Let us look at 
the spin dynamics at hi = to obtain information required to determine the parameters. The 
initial condition | ftt) leads (af (t)} = for any t. Let us take the initial state ^=(| t) + I 4 
))| |t)b instead, to avoid this problem. 
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Fig. A-l. (Color online) Eigenenergies as functions of hi when J\ — 4, J 2 = 5, h 2 — 6, and h 3 — 7. 



The Hamiltonian is block-diagonalized when h\ = as 

n u o 
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o n d 



J\ + J2 —hs -h 2 
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(A-9) 



-Ji + 


The Hamiltonian has four eigenvalues ±ei,±en, each of which is doubly degenerate. Here 

ei = y/s + 2y/C^, en = ^S -2^C~ (A-10) 

with 
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So = h 2 2 + h\ + Jl + J 2 2 , Co = fc|>i§ + /i§ Jl + J x 2 J 2 2 . 



(A-ll) 
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By using projection operators, time-evolution operators are expressed as 
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Then, real-time dynamics of x-component of the first spin is obtained as 
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